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ABSTRACT 


The major environmental Ultra Low Frequency (ULF), magnetic 
noise sources observed on the airborne AN/ASQ-81(V) optically 
pumped magnetometer are electromagnetic micropulsatrons, geo- 
logical structures and ocean wave noise. 

Experimental determinations of the amount o? magnetic noise 
produced by ocean surface gravity waves were made. The author's 
current research designed to remove the wave noise and increase 


the probability of Magnetic Anomaly Detection (MAD) is summarized. 
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Magnetometry has been used as a classical geophysical pro- 
specting tool for years (Dobrin, 1952). The principleg of 
Magnetic anamoly analysis has also been applied to Anti- 
Submarine Warfare (ASW) with a limited amount of success for 
submarine localization. 

The sensitivity of magnetometers has been continually 
improving so that it is now possible to observe, from an air- 


2 a 
nano Tesla/nertz Ἵπο ασ, imposea 


Bart, fluctuatious of 10° 
on the earth's constant field of approximately 50,000 nT. 

Land stationed superconducting magnetic snesors have demon- 
strated a sensitivity of Moe nlHz-phecddvenunt ob arbol 
superconducting magnetometers would greatly increase detection 
ranges* however, the extent of improvement depends in detail 
on the amplitude and frequency of the environmental noise which 
is above the noise floor of the instrument. 

Airborne Magnetic Anomaly Detection (MAD) for ASW is done 
in the .04 to 2 Hz frequency band. Ocean surface 
DIS the «03 to 10 Hz frequency band (Silvester, 1974) are con- 
tributors to the magnetic noise spectrum. The other two major 
environmental noise sources are electromagnetic puisations and 
geological noise. 

It is generally assumed that noise due to ocean gravity 


waves is not "loud" enough to be of any consequence. However, 


experimentally it was observed that noise levels induced by 





typical sea condition in the mid latitude oceans during the 
late fall to early spring are at least one order magnitude 
greater than the noise floor of airborne magnetometers. The 


Moise greatly reduces the effectiveness of Magnetic Anomaly 
| 


recognition. 
| 
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The problem of determining the electric current induced 


by deep water gravity waves at the ocean surface and the resul- 





tant periodic fluctuations in the magnetic field has been de- 


N 


eribed by Weaver [1965]. The magnetic field fluctuations 


above and below the ocean were obtained by solving Maxwell's 


- a 


Equations for the magnetic flux density vector in air and water, 
and applying electromagnetic boundary conditions at the sea 
Surface. 

The movement of the sea surface through the earth's mag- 
netic field induces electric eddy currents in the sea water. 
In deep water, assuming linear ocean wave theory, the water 
particle motion is periodic and circular, as is the induced 
magnetic field. This process becomes more complicated as one 
decomposes the sea surface spectrum into waves of various 
frequencies which may be propagating in various directions 
with respect to the magnetic field vector. Further complica- 
tions arises as higher order wave theories are applied to de- 
scribe the water particle motion and the waves propagate into 


Shallow water. 


ES MATHEMATICAL MODEL 
The mathematical model of the induced field was developed 
by Weaver [1965], assuming linear wave theory to describe the 


Sr particle motion. 





1t was derived for a monochromatic ocean surface gravity 
wave, propagating in deep water, in the X direction. The 
development is summarized below. 

The coordinate system formed a right hand triad as follows: 

X-axis, was taken in the direction of wave propagation. 

Y-axis, was taken parallel to the wave crest. 

Z-axis, was taken downward into the sea, Z-0 at the 
Still water level (SWL). 

It was assumed that the earth's magnetic field was constant 
in the region, and the wave crest was infinitely long. The 
model became two dimensional as all magnetic gradients in Y 


were negligible. 


. è — 
The earth's vector magnetic field, 3_, Was expressed as 


۸ 


A i‏ -- ^ ہت 
B, =[B.|(cosT cos6 £ - 9051 9.ν Sj+ sin L R)‏ 


^ ^ ^ 


Where 1, j, k were unit vectors along the three coordinate 
Beis x, y, Z. I, was the magnetic dip angle in the vertical 
plane, relative to the horizontal. 8, was the inclination 
in the horizontal plane of the direction of wave propagation, 


clockwise from the magnetic meridian (see Figure 1). 


The field vectors satisfied Maxwell's equation in the form 


2. 28 
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Volts 


کے کے 
where E was the electric field vector ( /Meter), B was the‏ 


induced magnetic flux density (Tesla), H was the permeability 


E iree space e ao Πε ο. O was the electrical con- 


Siemen 


کے 
Buctivity ( /Meter), V was the water particle velocity‏ 


Meter (Farads 


vector ( /sec), e was the permittivity /Meter), and 


t was time (seconds). 
The deep water particle velocity, described by linear wave 


theory, was of the form 


7 کے‎ ¿(ae -mx)- my Co 
V = ae er b.)e 


where a was the wave amplitude (Meters), w:27/f was the angular 


Radians 


frequency of the wave ( /sec), m-2m/L was the wave number 


E Trans eter), and t was time (seconds). 
—_ 
Since V was harmonic in x and t, the magnetic and electrical 


حم 
field resulting from V has a similar structure, thus‏ 


DS L(&E- mx) 
B= bia) 2 (4) 
جہ‎ = L(Wwt -mX) 


7 (5) 


B 


==> 
The B field and the field are non-divergent which im 


db EA 


rj 
E 
| 
H: 
(D 
A, 


= LM b.) (6) 
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EN L m Ext) ۱ (7) 
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The differential equation in water satisfied by B was 
obtained by assuming B««B, and we<<o. Then, taking the curl 


of equation (2) and using the non-divergent portion of the 


vector identity, Vx (VxA) = η + V(V-A), the following was 
obtained: 
- V^B -GYXEÉ-« SAMT X(VXBS . (8) 


Further simplification of equation (8) was obtained by the 
elimination of VxE by equation (i). From equation (3), V was 
eliminated. Then if y = Ίμοωσ, the following equation results 

2 950) ) L Y ( τ. η 

5 —+ = (y emi) Lum am ND oo - e mite Rg : (3) 
In region Z>0, equation (9) was solved for the Z component 


by assuming 


dr 
πη = Z D 2 4 ; ; — 
Data) = Q2 be / 2 am Ba- وروی‎ 2 


2 
where Q was an undetermined coefficient. 
In air £>>0, thus equation (2) became, 
= SE 
B . ςοξ 
YX = ES (10) 
ct 


Bae curl of euation (END when E was eliminated by equation 
(1), resulted in 
3 τι” 
à AC ος) 
-(m q En b ۱ (11) 


Confining the sea surface to scales of motion described by 


gravity waves implied that n^»»y, eu, thus 


123 


2 





ο ρω 


2A 


: m boy QUAD 


EG che region Z<0 the solution for the Z component of 


| equation (12) was of the form 
b (a) = m S 
b Re 


Ehere R was an undetermined coefficient. 


The magnetic field vector was continuous across the air-sea 


ct 


interface, thus, the following boundary condition at Z=0, existed: 


Dy az (2-0) - b, waren, (Z =O) 


CUL 7 


Ob, ae. (2=0) _ 2b; waren (Z=0) 


25 23 (14.5) 


Applying the above boundary values to equations (3) and (13) 


yielded 


R=Q+ am (Bai = L Be?) 


and 


2. \ 
7 کت‎ mit) z à. (B, 65 - L Bei) 1 57 


Equations (15) and (16) lead to 


Q= -2m a ( Reio -i δες») (17) 


fm t 6 y ^ wy 
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and 7 7 2 
R= am (Bein) ος Bay) mm) 
pe Cire) 
2 
Thus the Z component of the induced magnetic field vector 
| in alr was y Z 
: m 2, và کے‎ IM 
bat) = Am (8,00 e Boop Cx tm ) m) £ 
j ت‎ 
oe 
y 
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επ order to further simplify equation (19), 7<<l was true 





m 
for gravity waves resulting in 


C20) 


B c 
bj i i (Baco - ¿Beda "9 


From equation (6) and (20) it was seen that 
ch SD (021 
xy) Bix) 
Me weas convenient tOo measure the earth's field with a 


total field magnetometer. Thus, 


R te = +b 6A) Ces Tees tb s I . (22) 
MEASURED 

When equations (20) and (21) were used to simplify equation 
(22), and oniy the perturbations were considered, it can be 
shown that the magnitude of the observed magnetic field at a 


Specific frequency would be, 


a ہے‎ 


ἅπαν... 


b) - E IB. (ees oc (iz) ο ὔ ر‎ ον 


ZO, 


es 
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The deep water wave length (L) of a surface wave is related 


to the wave period (T) by, 
τ 
Ls al (24) 
ATT 
where g is the acceleration due to gravity. 


Eliminating m, from equation (23) by the use of equation (24) 


leads to, 


PEE os m 2, IN Απ. 
ح22۸ سا‎ (R (esT cose) sin) )p 37:2 = 


Equation (25) was used to model the magnetic field spectrum 
q S E 
when the directional ocean wave spectrum was known. 

Using an airborne magnetometer to measure the magnetic 


Spectra, there was a doppler type frequency shift, described as, 


p = #(i - (2mfu/g) es #) > 


c 


(see Appendix C). Where f' was the wave frequency observed on 
9 -1 = s. 
the aircraft (sec ^), I was the actual wave frequency, u was 


ne velocity of the aircraft Es 


/sec), and Y was the angle 
between aircraft motion and direction of wave propagation (see 


Figure 1). 
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Me TAE EXPERIMENT 


Magnetic data was collected by an airborne magnetometer 
flown over the sea surface in order to verify the theory pre- 
viously presented and to further understand the observed environ- 
ment in the 0 to 2Hz band. The AN/ASQ-81 magnetometer, aboard 
the P-3C Orion aircraft, is capable of measuring the signal 
produced by ocean gravity waves which are frequency shifted 
due to relative motion between the wave and the aircraft flying 
at 100 meters altitude. The analog voltage time series output 
Of the magnetometer was recorded on an HP3960, FM tape recorder 
for post flight analysis. 

The fluctuation of the sea surface below the aircraft was 
determined by two techniques. The analog voltage time series 
output of the APN-194 radar altimeter, also aboard the aircraft, 
was recorded. 

Independent of the aircraft data the ocean surface wave 
climate was determined from analysis of the Fleet Numerical 
Weather Central (FNWC), spectral Ocean Wave Model (SOWM) and 
measured wave spectrum from National Oceanic and Atmospheric 
Administration (NOAA), Environmental Buoy, NR-6003 (EB-17), 
located at 52.2N, 155.8W from the airborne data. 

The airborne data was collected in vicinity of SIE 156°W, 


E und the buoys’ reporting period, 00Z, 10 Feb 78. 


E AIRBORNE DATA COLLECTION 
The AN/ASQ-81 employs the atomic properties of optically 


pumped, metastable helium atoms to detect variations in total 


17 
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ME Ure 2. Region of data collection. (X) FNWC Grid Points, 
(A)NOAA Buoy, (---) First Leg of Aircraft Data, 
(-.-«) Second Leg of Aircraft Data. 
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field magnetic amplitudes. The output of the sensing element 
is converted to an analog voltage which is bandpassed and scaled 


Lip) as a funcion or time. 


to display field variations (Order 10 

Instrument specifications are described in Appendix A. 
Perpendicular legs, 270? True and 000? True, were flown at 

approximately 60 meters altitude and 180 knots (91.44 M/s) for 


30 minutes each. (See Figure 2) The magnetic data was collected 


En the .04 Hz to 2 Hz band, thus, eliminating the earth's con- 


Meant field and any micropulsations above 2 Hz. Aircraft motions 
were compensated for by the magnetometer using baroaltimetry, 
selective filtering and accelerometers. No variations in the 


magnetic field due to the aircraft were expected due to the 
above real time processing and the straight and level flight 
paths flown. 

The instantaneous aircraft altitudes were recorded simul- 
taneously with the magnetic data from the APN-19% radar altimeter, 
which responded linearly to apparent altitude changes due to 
long period surface gravity waves. The wave frequency was 
shifted due to aircraft motion (see equation (26)). Beam 
characteristics, S9 cone, were such that at 60 meters the 
diameter of the radar spot on the surface was approximately 70 
meters. ln order to resolve wave amplitudes from 60 meter 
altitude, the wave period must be greater than 10 seconds. 

The frequency and amplitude errors attributable to the size 
of the radar spot were considered in the analysis. 

The AN/ASQ-81 and APN-194 signal pick-off points are 
described in Appendix A. 


E) 





mo. SPECTRAL WAVE DATA COLLECTION 

The experiment was conducted around synoptic reporting 
period 002, 10 Feb 1978 Iimsvieınlty of NOAA Buoy, EB-17, so 
the spectral information computed by the Buoy for the 007 time 
f period was useful in interpreting the aircraft data. 

In addition to Buoy data the SOWM output was used in order 
to infer directionality and improve the spectral wave estimate 
Mi the region of the flight. 

The Buoy is designed to report meterological and oceano- 
graphic conditions at the sea surface every 3 hours. The wave 
conditions are measured for 20 minutes by accelerometers and 
reported for the Buoy location as frequency and amplitude 
(Steele and Johnson, 1977). The Buoy data of 00Z, 10 Feb 1973 


1s reported in Table I. 


Frequency (Hz) Period (sec) Amplitude (meters) 
ale qo Πεν 
07 ph vio dt 
21018 12 9 1.6064 
¿US ٤ 1.9505 
ai 100 deo +2 
ΠΕΙ 1 - ۷ 
a SS 02000007 
15 ΤΙ . +3 
¿14 T E 
15 2 . 6543 
“18 DS .Hh 793 
۷ھ‎ 5.5) eo del 
8 ت5‎ a La 

Table 1 


NOAA Data Buoy 00Z, 10 Feb 1978 
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The FNWC, SOWM uses the analyzed surface pressure field 
to determine wind velocities and then builds, radiates and dis- 
sipates waves as prescribed by wave theory (FNWC Tech Note No. 
ES). Model output is energy as a function of frequency 
and direction at each grid point. The wave energy (E) at a 
frequency can be converted to wave amplitude (a) by a= E-Af. 
SOWM frid coverage of the data area is shown in Figure 2. 


Bhe model output for the closest grid point is summarized in 


Eable II. 
DIRECTION FROM 29e autor 3u Os 
en) 
RIOD (SEC) AMPLITUDES (METERS) 
τε ES 0 کر‎ ΠΠ 
i .o3 REUS .u08 1.002 sue 
13.89 08 λος: 10039 ‚347 
12.82 eas 0 BE ος 
Mm? 04 0+000 ος Y .959 
10.99 της 1.584 ‚914 Tann 
2-71 1203 012 944 e 
8.5 944 eect ey τα; η 
TN ‚851 να ον, ση 
55 409 457 ach ‚167 
Table II 


ENWC, SOWM Data, 032 10Feb78. Grid point located at 
52.357N, 154.099W. Amplitudes (meters) 
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IV. DATA ANALYSIS 


η... ΓΞ SERIES ANALYSIS 

The general form of the magnetic and altitude time series 
is shown in Figure 3, which is a sample of each time series. 

The time series were digitized at .1 second intervals, 
resulting in a Nyquist frequency of 5 Hz. Any spectral power 
above 5 Hz would be folded back into the spectral estimate 
falliasing). The magnetic data, due to the low pass, 2 Hz, 
characteristics of the instrument, alliased only instrument 
noise which was at a low level. It was assumed that no sig- 
nificant alliasing of the altitude data occurred due to the 
altimeters 0.1 second RC constant which results in a 1.59 Hz 
putoff frequency T - Dro 

Prior to computing power spectral estimates, the data was 
calibrated to insure constancy of unit magnitudes. The magnetic 
data voltages were linearly calibrated to nT and the altimetry 
voltages were linearly calibrated to meters after the flight. 

The data was prepared for the Fast Fourier Transform (FFT), 
(Bendat and Piersol, 1971) by forming 16 ensembles, 102.4 seconds 
each. Each ensemble was then linearly detrended. The detrending 
removed significant variance from the frequencies below .001 Hz. 
The cosine function was used to taper 10% of the time series 
ey multiplying the first and last 51 points of the time series 
by the cosine. Points 1 and 1024 become zero and points 52 


thru 973 are as digitized. The tapering procedure is done in 


order to prevent non-real low frequency power leakage, due to 
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Horizontal 


Sample of magnetic time series (top) and 


altitude time series (bottom). 


۰ 5076+ط 
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ensemble averaging, into the spectrum. The variance removed 
in the time series taper is restored to the power spectra after 
ene FFT (Bendat and ۳ ٢٣٢٢٢ 1971). 

The FFT computes the real and imaginary components of the 
Fourier Transform. It was then possible to compute the sample 
power spectra of each signal by ensemble averaging. The normal- 
ized cross power spectra, defined as coherence and the phase 
of the cross power spectrum between signals were also computed. 

Each of the two flight courses were treated as separate 
data sets. Note in Figure 4 the power spectra of the altimeter 
(M*/Hz) and magnetometer (nT* /Hz) and the coherence and phase 
Bs a function of frequency. The numerical results are summarized 


in Appendix D. 


B. SPECTRAL WAVE ANALYSIS 

The analysis of the FNWC, SOWM directional wave spectra 
and EB-17 ocean wave amplitude spectra in order to describe 
the magnetic field variations observed aboard the aircraft, 
Bequired fitting the large scale information to the microscale. 
This change of scale was feasible due to spatial uniformity 
of the ocean wave spectra over the area and frequency range 
investigated. 

The ocean wave climate in an area can be described by that 
at a point over a period of time. EB-17 reported the waves 
during a 20 minute period to have the spectral shape listed in 
Table I. The FNWC ocean wave amplitude and directional spectra 


had a similar shape listed in Table II. 
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The directionality of wave components at EB-17 was deter- 
mined by group velocity analysis (Silvester, 1974) which pro- 
pagated the FNWC, SOWM data at nearby grid points to the buoy 
position. The greatest weight was on the grid point closest 

Ito the position of the buoy. The grid point was located at 
52.357N, 154.099W, 60 miles from the buoy. The analyzed buoy 


directional amplitude spectra is reported in Table III. 


PERIOD (SEC) AMPLITUDE (METERS) DIRECTION FROM (OT) 
16.7 .8322 34 09 
14.3 ‚7574 πρ 
12,5 1.6064 315° 
Π.Ι 1.9604 πο 
10.0 1.9312 320° 

ail 1.5419 220. 
8.3 1.0116 EU e 
7.7 .9083 SS 
Da ΠΟΙ; 330° 
6.7 66543 SN 
6.3 az 320° 
5.9 4511 SIS 
5.6 4212 So - 
Table III 


Analyzed directional wave data at NOAA buoy E3-17, 00% 10 Feb 78. 


It is important to note that the NOAA data and the FNWC 


L data have similar spectral shapes although amplitudes are not 


equal, verifying scale analysis assumptions. The buoy reports 
a Spectral peak at approximately 11 seconds.  FNWC predicts a 


Spectral peak at approximately 12 seconds. 
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جت 


The observed ocean wave climate on the aircraft, although 
similar, was not the climate observed by EB-17 or predicted 
by SOWM. Each of the ocean wave data sources had errors asso- 
ciated with them, which would contribute to the differences 
in the data. Scale considerations, length of airborne data 
Becord, modeled grid point locations, buoy location and flight 
path, further contribute to the uniqueness of the separate 
data sources. 

Due to the multidirectional nature of the SOWM data and 
its proximity to the flight path, it was assumed that the 
closest FNWC data point (Table II) most closely resembled the 
ocean wave climate observed on the aircraft. The SOWM wave 
information was applied to equation (25) resulting in an ex- 
pected magnetic amplitude spectrum. The expected frequency 
shift due to the aircrafts motion relative to the waves was 
Bomputed by applying equation (26). To be consistent with the 
power spectral density computed from the aircraft data the 
amplitude spectrum described above was squared and normalized 
by multiplying by the frequency interval of each frequency 
estimate. The results are summarized in Figure 5 and numeri- 


Bally in Appendix D. 
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Baker and Graefe (1968) showed experimentally that the 
theory reviewed in Section II describes the magnetic field 
induced by ocean surface gravity waves. Comparison of the 
aircraft magnetic power spectra, Figure 4, and the expected 
magnetic power spectra computed from SOWM data, Figure 5, 
further verifies the electromagnetic theory described above. 

Figure 6, which is the superpositioning of Figure 4a on 
Figure 5a and Figure 4b on- Figure 5b,illustrates that although 
the aircraft and wave data magnetic spectra are not the same 
they are similar in amplitude and frequency. The differences 
observed in the low frequency end can, in part, be attributed 
to the assumptions necessary to the formulation of each magnetic 
power spectrum and the assumption that each data source is 
stationary in time and space. 

The FNWC data, although close to the flight path, was not 
coincident with the aircraft track (see Figure 2). Furthermore, 
determining the barometeric surface pressure field over data 
Sparse areas, such as the Gulf of Alaska caused inaccuracies 
in the SOWM wave height, frequency and directional information, 
which could have contributed to some of the differences observed. 

When studying the magnetic spectrum computed from the air- 
craft data it cannot be assumed that the only source of noise 
Was the ocean gravity waves. Although no significant electro- 
magnetic micropulsations in the 0 to 2 Hz frequency band, were 


expected, they cannot be ignored. 
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Production of magnetic noise by the aircraft's motion and 
onboard equipage was considered. Ideally, most aircraft pro- 
duced noises were compensated for in the magnetometers real 
time processing scheme and the flight was straight and level 
on each leg which would further reduce noise. However, expe- 
rience shows that the noise sources described above are not 
completely removed from the spectral analysis when long data 
records are required. 

Note in Figure 6a that both the aircraft and the analyzed 
SOWM magnetic power spectra indicate a broad coincident frequency 
band. In Figure 6b it is seen that a change in the flight path 
Bauses the band width of coincidence to compress, shift to 
higher frequency,and two distinct maxima are observed in the 
aircraft spectrum. One of which is nearly coincident with the 
SOWM magnetic spectrum. This indicates that a significant mag- 
metic noise source ls external to the aircraft and ocean surface 
wave dependent. 

Most significant, 1s the information implied from the analysis 
of Figure 4 and Figure 6. It is seen (Figure 4) that the co- 
herence is high and the phase relationship is good between the 
magnetometer and the altimeter in a small frequency band (.2 to 
-/ Hz). The trend observed in the phase relationship over that 
band can be attributed to signal delays imposed by the electrical 
Miatvacteristics of each instrument. 

This highly correlated frequency band falls in the same 


band in which the ocean gravity wave noise is expected. The 





Mime ter specifications indicated that it was accurate, + 1 
Foot and has a linear response to about 1.5 Hz. Thus it can 
be concluded that the altimeter senses the sea surface flucuations. 

Note the large spectral peak below .1 hantz on the altı- 

C) 

meter power spectra of Figure h. The peak is approximately 
one order of magnitude greater in the data of the first flight 
leg (Figure 4a) when compared to the second (Figure 4b). This 
peak was due to an intermittent square wave caused by malfunc- 
tion of the altimeter. The square wave does not significantly 


degrade the analysis presented above. 
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VESSERPELICATION 


It was shown in section V that the altimeter can be used 
to describe the instantaneous fluctuations of the sea surface 
below the aircraft which are directly related to a portion of 
the magnetic noise. Thus, the high correlation between the 
altimeter and the magnetometer indicates that it is possible 
to remove a portion of the magnetic noise due to ocean surface 
gravity waves. 

Current Naval ASW tactics require that the magnetic field 
information be analyzed in the real time domain. Thus, noise 
cancelling must be done real time. 

Spectral analysis or Convolution Eiltering is noe practical 
due to thetsignal characteristics of the target anomaly and 
Ehe long length of records required to compute the FFT which 
describes the low frequency signals. 

The magnetic anomaly is a single transient signal imposed 


on the total magnetic field thus making spectral analysis 


Bmpossible. Furthermore, it is possible that the aircraft 


شس 


5070100 7 


fu 


May change direction and flight attitude before 
long record is collected from the altimeter and the magneto- 
meter to determine the coefficients of the Fourier Transform 
mer convolution filtering. 

In order to keep noise removal in the real time domain the 
principle of Adaptive Noise Cancelling (Widrow, et.al., 1975) 


Was applied to the 10 Feb. data. 
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The adaptive noise cancelling was done by a microprocessor 
developed at Naval Oceans System Center (NOSC) and already 
Ned in frequency analysis at frequencies higher than those 
analyzed for ultra low frequency magnetic information (0-2 Hz). 
The altimeter (reference signal) contained nolse that was corre- 
lated with noise in the magnetometer (primary signal). Each 
Signal was digitized and the algorithm used computed the solu- 
tion which removed the correlated noise from the primary by 
adaptively subtracting the reference Signal to converge on 
minimum power in the primary signal. 

Preliminary results indicate that this method may be useful 
for real time noise cancelling in airborne magnetic data (see 
Pigure 7 and Figure 8) while the magnetic anomalies not due 
to wave motion should be relatively unaffected. 

Note in Figure 7 the relative signal levels before and 
after noise cancelling. The signal appears to be cut in gain 
approximately four to one after noise cancelling. However, 
Spectral analysis (Figure 8) indicates that the only power 
reduction was in the frequency band that has been indentified 


as the ocean surface gravity wave noise frequency band. 
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Figure 8. Aircraft magnetic spectrum (A) compared to "noise 
cancelled" magnetic spectrum (B). When altimeter 
was used as noise reference. Each spectrum has 
16 degrees of freedom. Note the band limited 
Doise reduction indicated by the Gross πα σας 
area. 
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VII. SUMMARY 


It has been shown that ocean surface gravity waves con- 
tribute significantly to the environmental background noise 
observed in airborne magnetometry. The noise is the result 
of the magnetic field associated with electrical eddy currents 
induced by ocean wave propagation through the earth's magnetic 
preld. 

The APN-194 altimeter eboard the P-3C ORION aircraft is 
capable of describing the wave climate below the aircraft sim- 
ultaneously with total field magnetic measurements observed 
by the AN/ASQ-81(V) magnetometer, also aboard the P-3C. Thus, 
the altimeter can be used to infer the wave noise measured by 
the magnetometer. 

Analysis demonstrated that it is possible to remove a sig- 
nificant portion of the magnetico noise due to ocean gravity 
waves by using the altimeter as a noise reference and sub- 
tracting that processed reference in real time from the mag- 
netometer time series output. 

Further research with real submarine targets is underway 
in order that the ASW community will benefit from applications 
of this work. Development of à greatly improved magnetometry 


System is possible in the near future. 
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APPENDIX A 
MAGNETOMETER AND ALTEMETDTER SYSTEM CHARACTERISTICS 


me MAGNETOMETER SYSTEM PARAMETERS - 
AN/ASW-81(V) MAGNETOMETER GENERAL CHARACTERISTICS 


Performance-Magnetometer Group 


Range 23,000 to 75,000 gammas 
Noise level 0.016 gamma (peak-to-peak) 
Figure Of merit (FOM)? 0.24 gamma 


MAD bandpass (selectable) 


High Pass Conner requencies” O. 0006 95055 
and Verl Hz 
Rejection slope 36 dB/octave 
Low Pass Goruerrneggenoaxesc0.2909. (9x 
and 2.0 E 
Rejection slope 24 dB/octave 
Gain selection τ OU NR aml awe σος 208 
(full scale) and 40 gammas 
Altitude compensation Barometric type 


Reference altitude, self-zeroing 
Linear range, * 200 feet 


Environmental Limits 


This equipment meets the requirements of military specifica- 
tion MIL-E-5400 within the limits described below. 


Temperature 


(1) Equipment Operating: o E except the 
Magnetic Detector DT-323/AS5Q-81(V) or Towed Body 
TB-623/ASQ-81(V) which have +55°C maximum 
cem erature; 


(2) Equipment Nonoperating: -62?C to +85°C except the 
Magnetic Detector DT-323/ASQ-81(V) or Towed Body 
TB-623/ASQ-81(V) which have +71 C maximum 


temperature. 
Altitude Sea level to 30,000 feet 
0 1 O to 100 percent 
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SROCK 


Vibration 


Input Power Requirements 


| Magnetometer Group 


“ο το ο eg with duration of 11 
fee milliseconds 


Gta view) ον ο 6201500 Hz) “except 
Detecting Set Control C-6983/AS0-81(V) 
or C-9086/AS0-81(V) which are curve 
III (5 g, 5 to 500 Hz). 


115/200 volts ac 3-phase power 


Voltage limits 
Frequency limits 
Power 

Panel Lamp Power 
Voltage limit 
Frequency 
Power 


Mangetometer Pick off Points 


۰۷٠٠۶٠٦۶ 
300 το 20 ΠΩ 


150 VA maximum 


(Uto 295 VOS 
CE GOT a 


10 watts maximum 


enal piek oti posmnts xn che P-30. were Rack C-3, TB=433, 


pins C-14 and C-15. 
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mr. ALTIMETER SYSTEM PARAMETERS - 
AN/APN-194 RADAR ALTIMETER 


Manufactured by Honeywell, Inc., Government and Ared Products 
Division 


specifications 
μα ۰ 1 + +٤7۶٤ 


ES 


beam width to 3dB point 
.l sec RG constant 
25-35 nano second pulse 
0-25 nsec gate 
EE timeter Pick off Points 
Signal prek oft in the P-3C were FS 1005, T8-54]16, pins 


A-5 and A-6. 
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Figure B.-2 FNWC Northern Hemispheric 850mb Heights, 
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ARRENDEX C 
DERIVATTON OF SURFACE WAVE FREQUENCY SHIFT DUE 
TO MOTION OF AIRCRAFT RELATIVE TO THE WAVE 


The wave equation due to motion of the aircraft is, 


Lwt-mx~ iO te = Mee 
e = We 


Where w is the angular frequency of the wave, m is the 
Wave number, t is the time and x is the distance from a re- 
frence point in the direction of wave propagation. The primed 


terms are the values due to motion, 


a 


۶۳ یی +۶" 


Where u is the aircraft velocity ت3‎ is the angle clock- 
mee trom the horizontal wave motion vector to the horizontal 
craft motion vector. 


Then, 


,19t-m(x*(ucosU)t) | io t-mx 


Which leads to 


MI cos E t rt‏ 1ے 


then solving for w? leads to 


0^ = Wemu cosy 
when W = ٦ m = AT and L = μυς. 
equation 26 results 
ھ2‎ - = 27fu/ „cosy ] (26) 
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SPECTRAL VACUES OF AIRCRAFT DATA AND MODELED NOISE 
MAGNETIC CONSTANTS AT 52N, 156N FOR FEBRUARY 1978 
Earth's total magnetic field, Eur ۶ 7٣۲ 


Dip angle, T 


TO 





۳۰۳۷۰۴۰۰٠٠٠۹٣٠۶ ٦17۱۳۱۸ ت1٦‎ ۴۰۱۲۲ ۷۰۰۱۱۷ 78 ۰ 7 ۷ 
٥۳۱۶۰۰۲ ۲۰٠٦ 0.٥٣٢۳۰۲۳ "50۰۱۲۳ SPECTRA DATA LEG 1 


۶ ۶ ۹ ۵۸ AMPLE TUDE AT u) 
eR 0.0072 
0 2010494 
2 .002H 
.14 20037 
«19 00,29 
es .0041 
20 ‚0014 
22 “ο. οἱ 
ت2‎ .00064 
6 U05 
227 .0028 
g OIL 
22 .0039 
2 .0041 
237 007.6 
.42 ο 8 
45 . 0048 
.47 .0054 
36 0052 
nd 00165 
eo -0018 

I ¿00100 


Table D-l.a The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Aircraft Course 
ο. 
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nT 


FREQUENCY (HZ2) AMPLITUDE ( ΗΖ) 
11 .00087 
INS τση; 
alge PIDE 
Po +۷ 
5 ΠΤ. 
21 «0014 
«29 00S 
ο +۳ 
29 .0088 
ο “00003 
2ئ‎ .0087 
36 210059 
2:3 d +3 
.H2 sw OSs 
US 20/095 
با‎ 7 .0034 
5 00.2 
2906 20022 
58 0002 
7 00057 
89 ھ0‎ 
2 200105 

IS "ھ2‎ ٤ 


Table D-1.b The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Alrcarft Course 
μας. 
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Table D-2 FFT Numerical Output of Aircraft Data (Course 270°) 
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Figure D-1 Altitude Spectrum (CS 000%). 
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Table -3  FFT Numerical Output of Aircraft Data 
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